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We consider quantum wire arrays consisting of GaAs rods embedded in 
Ala;Gai_2;As and disposed in sites of a square or triangular lattice. The electronic 
and hole spectra around the conduction band bottom and the valence band top are 
examined versus geometry of the lattice formed by the rods, concentration of Al in 
the matrix material, and structural parameters including the filling fraction and the 
lattice constant. Our calculations use the envelope function and are based on the 
effective-mass approximation. We show that the electronic and hole spectra resulting 
from the periodicity of the heterostructure, depend on the factors considered and 
that the effect of lattice geometry varies substantially with lattice constant. For low 
lattice constant values the minigaps are significantly wider in the case of triangular 
lattice, while for high lattice constant values wider minigaps occur in the square 
lattice-based arrays. We analyse the consequences of our findings for the efficiency 
of solar cells based on quantum wire arrays. 

PACS numbers: 84.60.Jt, 73.21. Cd,17.40.+w 

I. INTRODUCTION 

Semiconductor periodic heterostructures have an additional, superstructure-related dis- 
crete translational symmetry, the period of which is much larger than that of the atomic 
lattice. This causes electronic and hole minibands, with minigaps between them, to form 
near the bottom of the conduction band and the top of the valence band, respectively, in 
the electronic/hole spectrum of the structure. In superlattices, heterostructures periodic in 
a one dimension (ID), the electronic/hole band splits at the boundary and in the centre 



of the first Brillouin zone to form minibands separated by minigaps of width equal to the 
difference between the conduction band bottom or the valence band top, respectively, in the 
component materials of the superlattice.^ Similar effects occur also in quantum wire arrays, 
semiconductor heterostructures periodic in two dimensions (2D), but the opening of mini- 
gaps and their width, determined in a more complex way by factors which include structural 
parameters in addition to material parameters, can only be studied numerically. The effects 
in question are also observed in other periodic composites, such as photonic, phononic or 
magnonic crystals, in the dispersion of electromagnetic, elastic or spin waves. -'^i^*^'^'^'^'^ 

Aroused many years ago and persisting ever since, the interest in periodic semiconductor 
heterostructures stems from the possibility of fabricating systems with the conduction or 
valence band split into a desired pattern of discrete minibands and minigaps^'^'^; this 
could be achieved through systems of semiconductor dots or rods of properly designed size 
and shape (cylinders, cubes, prisms, rings or antirings)^>i^ii^>i^. The discrete spectra of 
electronic conduction or valence bands significantly modify the electron or hole transport as 
well as the optical properties of the structure. Used in photovoltaic and thermoelectric cells, 
semiconductor heterostructures allow construction of devices converting lighir^*^'^'^'^ or 
temperature difference^'^'^'^ into electric current with efficiency higher than that allowed 
by homogeneous materials. Although the two effects differ in nature, the electronic or hole 
band structure plays an important role in both. The numerous theoretical studies on the 
subject use various approximations and methoda^'^'^Ii^SiiSEi^iSi of determining the energy 
spectrum of 2D and 3D heterostructures. Here we use the plane wave method with the 
effective-mass approximation, a technique successfully applied to studying the electronic 
states in ID, 2D and 3D heterostructures with quantum dots and wires of different shape 
and size, as well as interdiffusion and strain effects on electronic bands.— 

The main objective of this study is to investigate the effect of the crystallographic struc- 
ture on the electronic/hole spectrum in quantum wire arrays. In particular, we are going 
to determine the conditions conducive to the opening of minigaps in the conduction band 
and in the valence band. Besides its purely scientific aspect, our investigation is of practical 
importance for photovoltaic cells, the efficiency of which is determined by the electronic 
structure of the active material in which the photoelectric effect takes place. If the active 
material is a periodic semiconductor heterostructure, the efficiency of the cell will substan- 
tially depend on the width of the electronic energy minigap between the lowest miniband 



and the rest of the conduction band.— 

The paper is organized as follows. In Section [III we present in detail the plane wave 
method applied to the equation for the electronic envelope function in quantum wire arrays, 
and discuss the approximations used. Section [Till discusses the results of our calculations 
obtained for the conduction band in a quantum wire array consisting of GaAs rods em- 
bedded in AlGaAs and forming a square or triangular lattice. The spectra of these 2D 
heterostructures are studied versus Al concentration in the matrix material, filling fraction 
and lattice constant. In Section \W\ we analyse the energy spectrum of the valence band in 
the same heterostructures. The paper is summed up in Conclusions, in which we indicate 
the implications of our results for the efficiency of solar cells based on quantum wire arrays. 



II. PLANE WAVE METHOD IN EFFECTIVE-MASS APPROXIMATION 

Let us consider a quantum wire arrays consisting of infinitely long rods of material A 
forming a square or triangular lattice and embedded in a matrix of material B (Fig. 1). The 
rod axes are oriented along the z direction. The rods are assumed to represent potential 
wells for propagating electrons and holes. 

Here we shall only consider arrays with rods of circular cross section. The effect of the 
cross-sectional geometry of the rods has already been studied in our earlier paper—, in which 
we established that the efficiency of solar energy conversion in the photoelectric effect only 
slightly depends on the shape of the rods and that their circular cross section is optimum 
in this aspect. The filling fraction /, defined as the area ratio of the rod cross section to 
the unit cell of the heterostructure, will be a measure of the rod size. The maximum filling 
fraction value corresponds to the situation in which rods touch the unit cell limits. For a 
square or triangular lattice of cylindrical rods with cross-section radius p the filling fraction 
can be calculated from the following respective formulae: 

f = and f = , (1) 

where a is the distance between adjacent lattice sites. 

Our calculations of the electronic band structure of 2D semiconductor heterostructures 
will be performed in the effective-mass approximation^i^i^, which means that we shall 
restrict our attention to the bottom (and top) of the parabolic conduction (and parabolic 



valence) band in the electronic (hole) spectrum of the constituent semiconductors. Such 
approximations prove to be very useful in the electronic band calculations for semiconductor 
heterostructures with various component materials, as well as for single dots.— '^'^'^'^ We 
shall assume the interactions between the single conduction band and the two valence bands 
are negligible; this allows independent determination of the electronic miniband structure 
in the conduction band (Sections HTl and UTTl in the single band approximation) and the hole 
miniband structure in the valence band (Section llVp . 

Electronic states are described by the Ben-Daniel-Duke equation:— 



-« ( V-?T^V ) + Ec{v) 



^e(r)=m(r), (2) 



with effective mass m*(r) (isotropic in a homogeneous medium) and effective potential Ec{r) 
determining the position of the conduction band bottom; r denote the position vector in 
the plane of the periodicity; \E'e denotes the envelope wave function of conduction band 
electrons. The constant a = (2mee) ~ 3.80998 (where me and e denote the free 

electron mass and charge, respectively) allows to express the energy E in electronvolts (eV) 
and the coordinates x, y in angstroms A. The material parameters: the effective mass and 
the effective potential, vary in space with the periodicity of the heterostructure: 

m*(r + R) = m*(r), 

Ec{r + K) = Ec{r), (3) 

where R is a lattice vector (R = a(m^, iriy) in square and R = a{'mx + Tny/2, \^my/2) in 
triangular lattice, mj, and rriy being integers). 

The expansion of the electron envelope function in the plane-wave basis and the Fourier 
expansion of the material parameters can be written as: 

G 

Ecir) = Y.E^e^^\ 

G 

w{Y) = l/m*{Y) = ^M;^e^^■^ (4) 

G 

where 0^ are Fourier coefficients of the periodic factor of the envelope function (which has 
the same periodicity as the quantum wire arrays) and k is a wave vector from the first 



Brillouin zone. Fourier coefficients of the material parameters, and w^, can be found 
analytically from the formula: 

F"^ = T7 I f{r)e-"^'"dv, (5) 



where /(r) and denote, respectively, the periodic material parameter (the effective mass 
or the effective potential) and the corresponding Fourier coefficient for a plane wave of wave 
vector equal to vector G of the reciprocal lattice (G = {2Tr/a){nx,ny) for square lattice, 
G = {27r/a){nx, {2ny — nx)/^/S) for triangular lattice, Ux and Uy being integers); V denotes 
the area of the periodic heterostructure unit cell. The explicit form of the Fourier coefficients 
([5]) for rods of circular cross section is: 

Fa/ - Fb{1 - f) for G = 0, 
{Fa - FB)2f^ for G ^ 0, 
where functions Fa and Fb are the values of potential {Ec,a and Ec\b) or effective mass 
inverse {l/m*^ and 1/m*^) in the cylinders and in the matrix, respectively, and Ji is a Bessel 
function of the ffist kind. 

The substitution of expansions (jlj) in the Schrodinger equation ([2]) leads to the following 
system of equations representing an eigenvalue problem, the solution of which yields Fourier 
coefficients of the periodic factor of the envelope function and the electron energy E-r^'^ 

J2[a{G + k). (G' + k) w^'-^ + E^'-^] = E4>f. (7) 

G' 

When a finite number N of reciprocal lattice vectors is used in the Fourier series (jl]) also 
the system of equations ([7]) becomes finite. We have solved this system of equations by 
standard numerical procedures designed for solving symmetrical matrix eigenvalue problems, 
and tested the eigenvalues found for convergence, as necessary in procedures of this type. For 
all the structures considered in this study a satisfactory convergence of numerical solutions 
of Eq. (jll) for ffist few bands proves to be attained with the use of 169 reciprocal lattice 
vectors. 



III. THE EFFECT OF THE CRYSTALLOGRAPHIC STRUCTURE ON THE 

CONDUCTION BAND 



The results presented in this Section have been obtained for rods of gallium arsenide 
(GaAs) embedded in a matrix of gallium arsenide doped with aluminium (Al^Gai_^As, 



where x is the concentration of Al ions). In this case the depth of the potential wells 
is determined by the concentration of Al in the matrix: growing matrix Al concentration 
will increase the effective potential in the matrix, and thus deepen the potential wells felt 
by conduction band electrons. We have used the following empirical formulae for a linear 
extrapolation of the material parameter values in GaAs and AlAs to estimate their values 
in the AUGai_^As matrix: Ec = 0.944x and m* = 0.067 + O.OSSa;.-^'^ 

The effective-mass approximation is justified by the occurrence of a direct gap at point F 
of the atomic band structure for low concentrations of Al in both the rods and the matrix. To 
meet this condition of approximation applicability, the Al concentration in the matrix was 
only allowed to range from to 0.35. Also calculated in the envelope function approximation 
and by the plane wave method, the electronic band structure (the conduction band) of 
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471] is in agreement with the results 



square lattice-based heterostructures presented in 
obtained in this study from Eq. ([7]). 

Figure [2] (a) and (b) shows the energy structure of conduction band electrons versus Al 
concentration in a 2D heterostructure with GaAs rods disposed in sites of a triangular or 
square lattice. The calculations were performed at fixed lattice constant a = 120 A and 
filling fraction / = 0.3. Increasing x results in deepening potential wells, in which electrons 
from first minibands localise and the interaction between electrons in neighborhood wells 
weakens. As a consequence the miniband shrinks and increase their energy, and the minigap 
that separates it from the quasi-continuous conduction band above grows wider. There are 
also some differences between lattices: in the triangular lattice case the first miniband is 
seen to occur almost throughout the considered range of Al concentration x (from x = 0.01 
up), while in the square lattice case the Al concentration must be above a critical value 
X = 0.07 for the gap to open. Note from Al concentration x = 0.13 up the first minigap 
is wider in the square lattice case than in the triangular lattice case. The second minigap 
(between the third and the fourth miniband) only opens at x = 0.22 and x = 0.20, for the 
square and triangular lattice, respectively, and is wider in the square lattice-based structure. 

Let us see how the band structure varies with filling fraction. Figure [3] shows plots 
obtained for matrix Al concentration fixed at x = 0.35 and for lattice constant a = 120 A 
in the triangular and square lattice case. As expected, the minibands are seen to shift down 
as the filling fraction grows due to widening of the wells and to become wider because of 
shrinking of matrix (barrier) material^. The filling fraction dependence in both structures 



is very similar in terms of ranges of minigap occurrence as well as filling fraction values 
corresponding to the maximum width of the first and second gap, / ~ 0.18 and / ~ 0.39, 
respectively (see Fig. Hj). Note, however, that in the case of rods forming a square lattice 
both the first minigap and the second one are wider, their respective widths being 0.17 eV 
and 0.1 eV, against 0.14 eV and 0.08 eV in the triangular lattice case. These small differences 
between electronic bands in square and triangular lattices in dependence on x and / can be 
explained as follows. 

A major difference between a triangular lattice and a square lattice is the number of 
nearest neighbours; in a triangular lattice-based structure each well (rod) has 6 neighbours 
(at a distance a), against 4 neighbours (at a distance a) in the square lattice case. A larger 
number of nearest neighbours implies more interactions, if exist, and wider energy bands 
for the same parameter values, hence minibands in the triangular lattice case should be 
slightly wider than in a square lattice-based structure with the same values of x, f and a 
(especially those with high energy, e. g. compare Figs [2] (a) and (b) or Figs[3](a) and (b)). 
Obviously, this affects the minigap width, in which differences are seen as well. Fig. HI Let 
us get back to the definition of filling fraction, ([T]), to determine the radii of the cylindrical 
rods corresponding to equal filling fraction in both types of lattice at equal lattice constant 
values. The cylindrical rods are found to have larger diameter in the square lattice case 
(Psquare ~ l-0746ptrianguiar) 5 which implies relatively lower energies of the lowest minibands 
associated with electrons in wells, and consequently, increased isolation of these states, and 
narrower bands, also conducive to wider minigaps. However, it should be kept in mind that 
the spacing between rods decreases as their diameter grows, which can lead to increased 
overlapping of the wave functions of adjacent rods, and consequently, to wider minibands. 
The impact of these two competing effects can be expected to depend on the lattice constant. 

Let us compare plots depicting the band structure versus the lattice constant in the square 
and triangular lattice case (Fig. [5]). In this dependence the electronic band structures in the 
two cases considered differ the most. In the triangular lattice case the first minigap is seen to 
occur throughout the considered lattice constant range, i.e. from 25 A to 150 A, while in the 
square lattice case the gap only opens above a ^ 45 A. For the second gap the critical lattice 
constant values are a = 87 A and 92 A in the triangular and square lattice case, respectively. 
In Fig. inithe width of minigaps in the square and triangular lattice-based heterostructures is 
represented by the respective shapes, i.e. squares and triangles. Filled and outlined symbols 



refer to the first and second gap, respectively. For low lattice constant values a single gap 
is seen to occur in the triangular lattice case; the width of this gap decreases monotonously 
with increasing lattice constant from 0.3 eV for a = 25 A to 0.06 eV for a = 250 A. In the 
square lattice case the width of the first gap is seen to vary with a in a non-monotonous 
manner: from a = 42.5 A up the gap widens to reach a maximum width of 0.16 eV at lattice 
constant 100 A, and then narrows down as a increases still further. For a = 75 A the first 
gap has equal width in both cases. The second gap has a maximum width for a ^ 155 A 
for both lattice types considered, but tends to be wider in the square lattice case. Thus, for 
large lattice constant values we observe the relations already explained above. 

Let us try to elucidate the intriguing relations found in the range of low lattice constant 
values. Two effects occur as a result of reducing the lattice constant at fixed filling fraction: 
(a) the rod (potential well) diameter decreases, with a concurrent increase energy of bands; 
at the same time, (b) the distance between adjacent wells decreases, which implies increased 
overlapping of wave functions and thus steps up the widening of bands. But it doesn't explain 
relations for a < 55 A where the top of the first miniband is already above the barrier, for 
both arrangements and electron waves can propagate across the structure slightly scattered 
by the barrier. It means, that triangular lattice, due to a larger number of adjacent wells, 
favors destructive interference of such waves resulting in wide gap, as opposed to square 
lattice, where the gap disappear. Figure [7] shows the band structures, calculated along the 
path r = (0, 0) X = 27r/a(0.5, 0) ^ J = 27r/a(0.5, 1/(2^3)) ^ F and F = (0, 0) ^ X = 
27r/a(0.5,0) — M = 7r/a(0.5, 0.5) — F in the first Brillouin zone, of triangular (solid lines) 
and square (dashed line) lattice-based 2D heterostructure, respectively, with lattice constant 
a = 50 A (a) and a = 120 A (b). As expected the bands are seen to be very wide in the 
quantum wire arrays with low lattice constant value. Fig. [7] (a), and very narrow in the arrays 
with large lattice constant Fig. [7] (b). In Fig. [7| (a) we can observe additional difference in 
band structures for triangular and square lattices. For triangular lattice separating the first 
miniband form the second one at X point remains almost unaltered in comparison to point 
J. In the case of square lattice the separations both in X and M point are reduced as lattice 
constant decreases. This effect and some small shift of miniband position between points X 
and M cause the closing of absolute minigap in square lattice whereas for triangular lattice 
this minigap remains open. 



IV. THE EFFECT OF THE CRYSTALLOGRAPHIC STRUCTURE ON THE 

VALENCE BAND 



For calculating the valence band structure of quantum wire arrays we use the plane wave 
method, applied to the Schrodinger equation of the envelope function of light and heavy 
hole states near the top of the valence band:— 




= E 




(8) 



where "^ihir) and ^E^h/ilr) are envelope functions for light and heavy holes, respectively, and 
operators P, Q and R have the form: 
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Luttinger parameters 71, 72, 73, describing the effective masses 1/(71 + 72) and 1/(71 —72) of 
light and heavy holes, respectively, near point F of the atomic lattice, are, like the position 
of the valence band top Ey, periodic in the heterostructure. 

Expanded in the plane-wave basis by a procedure similar to that employed in the case 
of the conduction band, equation (IHl) can be written as a Hermitian eigenvalue problem for 
interacting light and heavy holes: 



J2[iA, + A,) 0,^' + (5i - tB2) 



J2\{B,+ iB^) + (Ai - A 



G' 

2jVhh 



(10) 



where Ai, A2 and Bi, B2 are expressed as follows: 

A, = -a^^-""' (G + k) ■ (C + k) - E^-^', 
A2 = -0^2'""' (G + k) ■ (G' + k) , 

-{Gy + k,){G'^+ky)], 

B2 = aV?>if-'''[{G', + K){Gy + ky) 

+ + + (11) 

0^ and (f)^ have the sense of Fourier coefficients for periodic factors of the envelope functions 
of heavy and light holes, respectively, and 7p = 1, 2 or 3) and Ey are the respective 
Fourier coefficients of Luttinger parameters 7^ and Ey, which can be found from ^ for 
the Luttinger parameter values, specified below, in the rod material A and in the matrix 
material B. 

Again, let us consider a heterostructure consisting of a triangular or square lattice-based 
system of GaAs rods embedded in Alj,Gai_j.As. The following material parameter values, 
dependent on the concentration of Al in aluminium gallium arsenide, can be assumed:—!^ 

Ev = 1.519 + 0.75X, 

71 = 6.85 - 3.40a;, 

72 = 2.10- 1.42a;, 

73 = 2.90 - 1.61a;. (12) 

Fig. [8] shows the wave vector dependence of hole minibands plotted along the irreducible 
part of the first Brillouin zone in a triangular lattice-based heterostructure (solid line, path 
F X ^ J ^ r) and in a square lattice-based heterostructure (dashed line, path F — »• 
X ^ M F). The calculations were performed for lattice constant a = 50 A, filling fraction 
/ = 0.3 and matrix Al concentration x = 0.35. 

Comparing the hole energy structure (Fig. [8]) with the energy structure of conduction 
band electrons (Fig. [7j), we notice the hole bands (for light and heavy holes) are much 
narrower and arranged very densely, so their role in photoelectric effect will be minor. The 
first three electronic bands lie in the energy range from 0.21 eV to 1.14 eV, against -1.76 eV to 
-1.65 eV in the case of hole bands (the negative energy values are a consequence of assuming 



the energy Ec^a of a conduction band bottom electron as the zero energy level). The lattice 
type-related differences in energy structure will be better visualised in the dependence of 
hole minibands on the lattice constant and the filling fraction. 

The lattice constant dependence of the band structure in the triangular lattice case is 
depicted in Fig. [^a). The minibands are seen to shift towards higher energies and shrink 
as the lattice constant increases. A similar tendency is seen in the square lattice case, in 
which, however, the minibands are shifted up: for a = 250 A the shift is 0.03 eV. Narrow 
bands are characteristic of weakly interacting states localised in potential wells, and their 
large number is a consequence of the low valence band maximum {Ey^B = —1.78 eV) in the 
matrix material. The minigap width, plotted versus lattice constant in Fig. Mjo), is seen 
to grow steeply at first to reach a maximum and then reduce monotonously as the lattice 
constant continues to increase. This behaviour is typical of all the minigaps found in both 
the triangular lattice-based array and the square lattice-based structure. Minor differences 
in minigap position and width only occur for low lattice constant values and disappear as 
the lattice constant grows. 

Let us scrutinize the evolution of the hole minibands with filling fraction. Figure [TOT a) 
shows the highest two hole minibands versus filling fraction in the triangular and square 
lattice case (lines with triangles and squares, respectively). The bands are seen to shift 
towards higher energies with increasing /. In the square lattice case this increase in energy 
is much stronger, which implies wider minigaps, as depicted in Fig. [TOT b). showing the 
width of the two minigaps plotted versus filling fraction. The minigaps are the widest {AE 
up to 0.05 eV) for low filling fraction values, / ~ 0.05, and their maximum width represents 
less than one third of that of electronic minigaps in the valence band. 

Increasing frequency of holes minibands with increase of lattice constant and filling frac- 
tion is due to the same mechanism as has been already explained in Sec. Illll for the frequency 
of electron minibands decrease. This is because cylinders of GaAs forms wells for the holes 
in valence band, as well as for electrons in conduction band. The observed differences in 
both dependencies (for electrons and holes) are mainly due to the number of states localized 
in these wells. We can compare this numbers with the help of thick (red) lines in Figs [7] and 
in Fig. [H] for electrons and holes respectively. In the case of holes, the Alo.35Gao.65As form 
the barier with the height -1.78 eV and the large number of hole minibands have energy 
inside the wells. This is responsible for similar behavior of hole minibands and minigaps in 



the presented in Figs M and [TU] dependencies on lattice constant and filling fraction. 



V. CONCLUSIONS 

We have examined thoroughly the effect of the crystallographic structure on the elec- 
tronic and hole energy bands in quantum wire arrays by comparing structures consisting of 
GaAs rods embedded in Ala;Gai_a;As and disposed in sites of square or triangular lattice. 
We have established that both the electronic and hole minigaps tend to be wider in the 
square lattice-based structure in the considered ranges of parameters that include the lat- 
tice constant, the filling fraction and the potential well depth. Exceptions from this rule are 
heterostructures with low (less than 75 A) lattice constant values (with the first electronic 
minigap, between minibands 1 and 2, wider in the triangular lattice case than in the square 
lattice case) and heterostructures with low Al concentrations in the A^Gai-^^As matrix. 
We have demonstrated that the conduction band bottom and the valence band top can be 
controlled in the range ~ 0.3 eV and ~ 0.2 eV, respectively, by adjusting the filling fraction, 
and in the range ~ 0.14 eV and ~ 0.1 eV, respectively, by adjusting the lattice constant. 
Also the minigap width is shown to be controllable, especially the width of the gap between 
minibands 1 and 2, which can attain (at matrix Al concentration 0.35) 0.19 eV and 0.06 eV 
in the conduction band and in the valence band, respectively. 

The above-discussed properties of semiconductor heterostructures will be crucial for the 
efficiency of solar cells based on these materials. In monolithic semiconductors with a sin- 
gle bandgap the absorption of photons is incomplete because photons of energy below the 
bandgap width are not absorbed by the system. On the other hand, thermalisation causes 
the surplus energy of electrons above the bottom of the conduction band or holes below 
the top of the valence band to dissipate in thermal contact with the crystal lattice. Photon 
absorption can be increased by using tandem cells - cascades of solar cells of successively 
narrowing bandgaps^"^. A similar effect can be obtained with solar-energy converters based 
on semiconductor heterostructures discussed in this paper with multiple bandgaps, which 
also allow to reduce thermalisation losses.— 

AlGaAs heterostructures have relatively wide energetic gap (1.4-2.1eV) which make them 
potentially suitable for designing of solar cells with intermediate band, were the wider gap 
(e.g. in comparison to the bulk Si) is required for achievement of high efficiency of solar 



radiation conversion. The key role in this gain in efficiency is played by the lowest con- 
duction miniband, which, detached from the overlapping conduction minibands, acts as an 
intermediate band that opens an extra channel for carrier transitions between the valence 
band and the conduction band. Another parameter of vital importance for the efficiency 
of solar-energy conversion is the distance between the top of the highest valence miniband 
and the bottom of the overlapping conduction minibands. This distance, determines the 
energy of utilized carriers.— The effect of the other electronic minigaps on the efficiency 
of solar-energy conversion remains to be investigated. Considering that for ffiling fraction 
/ > 0.5 and/or sufficiently large lattice constant values (a > 160 A for / = 0.3) the sec- 
ond minigap in the conduction band is comparable to or even wider than the ffist one, it 
can be expected to affect the efficiency of solar-energy conversion by providing additional 
allowed transitions. However, a precise evaluation of the effect of further discretisation of 
the conduction or valence band is beyond the scope of this study and will be the object of 
our future research. The impact of the discretisation of the valence band on the efficiency of 
the photovoltaic effect will be much lesser than that of the discretisation of the conduction 
band, and, in the ffist approximation, can be regarded as limited to a shift in valence band 
top in the evaluation of solar cell efficiency.— 
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FIG. 1: Cross section of a 2D heterostructure consisting of infinitely long cylindrical rods (A) 
disposed in sites of a (a) triangular or (b) square lattice and embedded in matrix material (B). 
Dotted lines represent unit cell limits, a is the lattice constant, and p the radius of the circular 
cross section of the rods; in both depicted structures the filling fraction is / = 0.3. 
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FIG. 2: Electronic band structure of 2D superlattice formed by periodic system of cylindrical GaAs 
rods disposed in sites of a triangular (a) or square (b) lattice and embedded in Ala;Gai_2,As, plotted 
versus matrix Al concentration x at fixed filling fraction / = 0.3 and lattice constant a = 120 A. 
The thick (red) line indicates the conduction band bottom energy in the Al2:Gai_3;As matrix. 
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FIG. 3: Electronic band structure of 2D superlattice formed by periodic system of cylindrical GaAs 
rods disposed in sites of a triangular (a) or square (b) lattice and embedded in Alo.ssGao.esAs, 
plotted versus filling fraction / at fixed lattice constant a = 120 A. The thick (red) horizontal line 
indicates the conduction band bottom energy, Ec,b ~ 0.33, in the matrix. 



First gap, square lattice 




0.2 0.4 0.6 

Filling fraction, f 

FIG. 4: Width of first and second energy minigap (solid and outlined symbols, respectively) versus 
filling fraction / in a triangular or square lattice-based heterostructure (triangles and squares, 
respectively). Heterostructure form quantum wire arrays of cylindrical GaAs rods disposed in sites 
of a triangular or square lattice and embedded in Alg.ssGao.esAs; lattice constant a = 120 A. 
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FIG. 5: Electronic band structure of quantum wire arrays formed by system of cylindrical GaAs 
rods disposed in sites of a triangular (a) or square (b) lattice and embedded in Alo.35Gao.65As, 
plotted versus lattice constant a at fixed filling fraction / = 0.3. The thick (gray) horizontal line 
indicates the conduction band bottom energy in the matrix material, Ec\b- 
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FIG. 6: Width of first and second electronic minigap (filled and outlined symbols, respectively) 
versus lattice constant a in triangular or square lattice-based heterostructure (triangles and squares, 
respectively) at fixed filling fraction / = 0.3. The heterostructure consists of a periodic system of 
cylindrical GaAs rods embedded in Alo.35Gao.65As. 
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FIG. 7: Electronic minibands in system of GaAs rods disposed in sites of triangular (solid lines) 
and square lattice (dashed lines) with lattice constant a = 50 A (a) and a = 120 A (b). Rods are 
embedded in Alo.35Gao.65As, and the filling fraction is / = 0.3. Note the different energy scale in 
the two plots. The thick (gray) horizontal line indicates the conduction band bottom energy in the 
matrix material, Ec,b- 
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FIG. 8: Hole minibands in heterostructure with GaAs rods disposed in sites of triangular or 
square lattice (solid and dashed lines, respectively), calculated along the indicated path in the first 
Brillouin zone; bottom and top scales refer to triangular and square lattice-based heterostructure, 
respectively. In both cases the matrix material is Alo.35Gao.65As, the lattice constant a = 50 A, 
and the filling fraction / = 0.3. The thick (red) horizontal line indicates the valence band top 
energy, Ey^B = —1.78, in the matrix material. 
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FIG. 9: (a) Hole minibands versus lattice constant in heterostructure with GaAs rods disposed 
in triangular lattice sites, (b) Width of the first two hole minigaps, between minibands 1 and 2 
(filled symbols) and between minibands 3 and 4 (outlined symbols), versus a for rods disposed 
in triangular or square lattice sites (triangles and squares, respectively). The matrix material is 
Alo.35Gao.65As and the filling fraction / = 0.3. The valence band top energy is Ey^j^ — — 1.78 in 
this case. 
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FIG. 10: (a) Energy of the first two hole minibands versus filling fraction in quantum wire array 
with GaAs rods disposed in triangular or square lattice sites (triangles and squares, respectively), 
(b) Width of top two minigaps, between minibands 1 and 2 (filled symbols) and between minibands 
3 and 4 (outlined symbols) versus / in the triangular and square lattice case (triangles and squares, 
respectively). The matrix material is Alo.35Gao.65 As, the lattice constant a = 50 A, and the filling 
fraction / = 0.3. 



